Spontaneous helicity of a polymer with side-loops confined to a cylinder 
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Inspired by recent experiments on the spatial organization of bacterial chromosomes, we consider 
a type of "bottle brush" polymer consisting of a flexible backbone chain, to which flexible side loops 
are connected. We show that such a model with an open linear backbone spontaneously adopts a 
helical structure with a well-defined pitch when confined to small cylindrical volume. This helicity 
persists over a range of sizes and aspect-ratios of the cylinder, provided the packing fraction of the 
chain is suitably large. We analyze this results in terms of the interplay between the effective stiffness 
and actual intra-chain packing effects caused by the side-loops in response to the confinement. For 
the case of a circular backbone, mimicking e.g. the E. coli chromosome, the polymer adopts a 
linearized configuration of two parallel helices connected at the cylinder poles. 

PACS numbers: 82.35.Pq, 87.16. Sr, 87.15. ap, 89.75.Fb 



That confinement can dramatically influence the prop- 
erties of non- ideal polymers has already been appreciated 
theoretically for a long time. Scaling arguments suggest 
that when a confining dimension D becomes smaller than 
the natural size of a linear polymer, typically given by its 
radius of gyration R, the polymer, rather than behaving 
like a single coherent 'blob' of size R, effectively becomes 
a string of 'blobs' of size D [![, signalling a cross-over 
to lower-dimensional behavior. The experimental explo- 
ration of this regime, however, is difficult using synthetic 
polymers, whose typical radii of gyration are in the nm 
regime and moreover tend to be highly polydisperse in 
length 0. Fortunately nature provides us with an im- 




FIG. 1: (Color online) Helical equilibrium structure for a 
backbone chain of length h = 200a [black (blue) thick line] to 
which side-loops of length l s = 40a are attached at each back- 
bone monomer. The polymer is confined within a cylinder of 
length L — 50.75cr and diameter D = 29.5a. The side-loop 
monomers are shown as transparent gray (green) beads. 



portant class of model polymers in the form of DNA. To 



start with DNA is produced with almost perfect length 
control. A well-studied example in vitro is purified A- 
phagc DNA which contains exactly 48502 basepairs and 
has a length of 16.5 {im and a radius of gyration of 2.1 
/im. Using biochemical tools this length can be scaled 
to multiples of this unit. Moreover, DNA can readily be 
fluorescently labeled. These two properties of DNA were 
exploited in a number of recent experiments aimed at 
studying the behavior of polymers confined to nanoflu- 
idic channels 

More importantly, in vivo DNA in the form of chro- 
mosomes is almost invariably strongly confined. The 
paradigmatic case is the 1.5 mm long circular chromo- 
some of the bacterium E. coli, which is confined to a 
cylindrical volume of diameter ~ 0.8^m and length vary- 
ing between 2 — 4/im. Over the past few years there has 
been increasing interest in understanding the details of 
the chromosomal configuration in E. coli and its impli- 
cations for the, as yet far from fully understood, mech- 
anism of sister-chromosome segregation prior to division 
@~@|- Strikingly, recent microscopic observations have 
shown that bacterial chromosomes can exhibit a helical 
spatial density distribution within the cellular confine- 
ment, with a pitch-length in the order of a fraction of 
the cell length 0,[T3|. 

Here our aim is to explore to what extent the physics of 
highly confined polymers by itself provides clues to such 
novel global configurations. As a first step in that direc- 
tion it is necessary to take into account that the structure 
of a typical bacterial chromosome is not simply that of 
a closed linear polymer chain. The combined effects of 
supercoiling, due to a globally maintained under-twist 
of the DNA, the action of a number of chromosome re- 
modelling proteins, or even electrostatic "zippering" can 
cause loops in the DNA [ll[ , indicating that a more com- 
plex structural picture is needed. Although in reality 
the dynamics of the formation and the statistics of such 
loops are likely to be highly complex, almost certainly 
involving polydispersity in loop sizes and topological en- 
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tanglcments, we chose to focus on arguably the simplest 
model of a polymer "dressed" by a cloud of loops, i.e. a 
backbone chain to which side-loops of equal size are at- 
tached at a regular spacing. The same model has recently 
also been discussed by Reiss et al. [l2[, wno have stud- 
ied its behavior in the absence of confinement effects. 
This type of polymer model is similar to the so-called 
"bottle-brush" polymers, extensively studied by Binder 
and co-workers 13 - Il5j . The latter work has shown that 
such polymers develop a local resistance to bending due 
to the entropic repulsion between the side chains. This 
effective stiffening, combined with intra-chain packing ef- 
fects within the cylindrical confinement, leads, as we will 
show in the following, to the spontaneous formation of 
helical configurations. 

Our model chromosomes are of the bead-spring type, 
with consecutive beads attached to each other by a 
harmonic spring Vb = (A/2)(di — em,) 2 where di 
r, . i n, n is the position of z-th bead, a the equi- 
librium bond- length and Uj = di/|di| is the local tan- 
gent vector to the chain. Non-bonded beads repel each- 
other through the Weeks-Chandler- Andersen (WCA) po- 



tential [16j V WC A = 4e [{a/nj) 12 - (cr/r y ) 6 + 1/4] if the 



inter-monomer separation rjj < 2 1 / 6 er else Vwca = 0, 
where e and a set the energy and length scale of the 
system respectively. We use A — lOOe. The interac- 
tion of all beads with the confining walls are modelled 
through V waU = 2ire[(2/5){a/r lw ) w - ( ff /r„) 4 + 3/5] if 
the distance of the i-th monomer from a wall Ti W < a 
and V wa u = otherwise. We simulate this system em- 
ploying a velocity- Verlet molecular dynamics scheme in 
presence of a Langevin thermostat fixing the temperature 
at UbT = 1 as implemented by the ESPResSo package 

ill- 

We first consider a polymer composed of a linear back- 
bone chain of length 4 = 200a to which side-loops of 
length l s = 40cr are attached at every backbone monomer 
of the main chain. This polymer is confined to a cylinder 
of length L = 50.75ct and diameter D = 29.5c, yielding 
a monomer packing fraction of r\ — 23.8%. In Fig. [T] we 
show a typical equilibrium configuration of this polymer, 
which evidently displays a marked helical ordering of the 
backbone chain. The degree of helical ordering can be 
quantified by considering the tangent-tangent correlation 
function (u(s) • u(0)), where the positional coordinate is 
given by s = ia with i = 0, 1, . . . , 200. The Fourier trans- 
form of this quantity yields a structure function S(q) with 
a peak at a dimensionless wavenumber q max = h/^max, 
where X ma x is the pitch of the helix measured along the 
backbone chain. The height of the structure function at 
its maximum is a relative measure of the degree of he- 
licity, whilst the width of the peak is indicative of the 
statistical dispersion of the structure. 

Fig. [2] shows the correlation function and the corre- 
sponding structure function for this polymer, as we vary 
the diameter of the confining cylinder. This shows that 
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FIG. 2: (Color online) The tangent-tangent correlation (inset) 
and its Fourier transform for the polymer shown in Fig. [T] 
The correlation function is oscillatory, with the periodicity 
captured by the peak in the structure factor at q p = 4 for 
D/a = 14.68, 19.85, 29.5, and q p = 2 for D/a = 58.93. Also 
shown are the corresponding results for a main-chain polymer 
with an effective Gaussian-core mimicking the effect of inter- 
side-chain repulsion (cf. Fig. HI&)). 



the helical pitch is relatively robust against changes in 
the diameter, although a slight decrease in the amplitude 
is apparent as we increase the diameter. Only for the 
largest diameter, when both the degree of confinement 
as well as the overall packing fraction are significantly 
decreased, do we see a preference for a more longitudinal 
packing of the main chain. 

We now argue that the helical arrangement is stabilized 
by two effects. The first is the effective stiffness induced 
in the backbone by the presence of side loops. To quan- 
tify this intrinsic effect in the absence of confinement we 
study a free polymer of length lb = 500, with side-loops 
of the same length as before, l s = 40a, grafted at each 
backbone monomer. Although for a very long backbone 
(lb 3> l s ) one expects the exponent b, which governs the 
mean-square separation between monomers at a distance 
s along the backbone through the scaling (r(,s) 2 ) ~ s 2b , 
to reproduce the value b = 3/5 of a simple self- avoiding 
polymer, we find at the shorter length-scales relevant to 
our confined polymer a much higher value of b = 0.97 
(FiglH inset). At the same time, the tangent-tangent 
correlation (u(s).u(O)) (FigJ3|) shows an extremely weak 
power-law decay s~ a with a = 0.06, which satisfies the 
generic scaling rule a = 2 — 2b. The small value of the 
exponent a <g; 1 suggests that at these relatively short 
length scales the backbone is characterized by a signifi- 
cant effective stiffness, even close to that of a of a rigid 
rod for which 6=1. The algebraic nature of the corre- 
lation decay, however, precludes a simple interpretation 
in terms of an intrinsic length, in contrast to the persis- 
tence length of a chain with intrinsic resistance to bend- 
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FIG. 3: (Color online) The tangent-tangent correlation shows 
a weak power-law decay s~ a with a = 0.06. Inset: power-law 
growth of distance between two segments R(s) 
s b with b = 0.97. 



ing. Nevertheless, it is intuitively clear that there is a 
free-energy cost associated with local deformations of the 
backbone caused by the changes induced in the side-loop 
packing. For comparison, we note that the end-to-end 
distance measured for our backbone with If, — 500, which 
we determined to be R « 413 a, would be reproduced by 



a worm- like chain with persistence length P pa 391 a [18 1 
comparable to the contour length. 

However, backbone stiffness alone is not sufficient to 
explain the emergent helicity. Although earlier work has 
shown that a persistent chain confined to the surface of a 
cylinder can adopt helical confirmation [l9|], a persistent 
chain confined to the volume of a cylinder does not. We 
confirm this by simulating a worm-likc-chain (WLC) of 
length If, = 200cr with a persistence length of P = 2Z5, 
close to the effective value determined above for the free 
polymer. A typical equilibrium structure of the confined 
WLC is shown in Fig. 21 which clearly displays the ten- 
dency of such a chain to align itself with the long-axis 
of the cylinder, without much internal structure develop- 
ing, also supported by an analysis of the tangent-tangent 
correlations (see supplementary material). 

This negative result for the persistent chain, points 
to the importance of the packing of the side-loops in 
stabilizing the helical structure. One could naively ex- 
pect that the interactions between the side loops cause 
an additional effective soft repulsion between the back- 
bone monomers with a range determined by the radius of 
gyration of the side-loops. This repulsion has a dual role; 
as we have shown above, it stiffens the backbone at short 
length scales, but at long length scales it acts to keep 
distant parts of the backbone apart, effectively thicken- 
ing the backbone to a soft 'tube'. To assess the validity 
of this latter idea we consider a chain whose monomers, 
apart from interacting through the short range WCA po- 





FIG. 4: (Color online) (a) Self-avoiding WLC and (6) self- 
avoiding Gaussian-core polymer within a cylinder. The chain 
has contour length li, = 200cr, and for the WLC polymer the 
persistence length P = 2lt- The confining cylinder is of length 
L — 50.75cr and diameter D = 29.5(7 for both (a) and (b). 



tcntial, also interact with an effective Gaussian core in- 
teraction 



0V gc (r) =aexp[-(r/ W ) 5 



(1) 



intended to mimic, as simply as possible, the soft repul- 
sion between the side-loops. Bolhuis et al. showed that, 
in free space, the effective interaction between two linear 
polymers of the same length can be approximated by the 
above expression with a ~ 2 and w = R g , the radius of 
gyration of the individual polymers [2^ . While it is well 
known that the free energy cost of overlap between two 
polymers in free space is independent of polymer size [EL 
I23J , the blob picture of polymers suggests that their over- 
lap free energy will become a linear function of polymer 
length Q when the chains are strongly confined. We 
therefore take the interaction strength a = 40 to mimic 
the densely-packed side-loops of size l s — 40er. From our 
simulations of side-loop coupled back-bone, we found the 
mean radius of gyration R g = 2.95 a of the side-loops of 
length l s = 40cr. However, due to the repulsive inter- 
actions with the backbone monomers the center of mass 
of the loops is offset from the backbone. We therefore 
extracted the distribution of the center of mass of the 
side-loops with respect to their attachment point on the 
backbone in the plane perpendicular to the local tangent 
direction along the chain (see Figure [5] in the Supplemen- 
tary Material). This distribution has a sharp maximum 
at a distance of approximately 4er from the backbone 
chain. This result also shows that indeed the density dis- 
tribution of the side-loop monomers is structured in a 
manner consistent with the 'thickened tube' picture al- 
luded to above. We therefore took the size parameter 
characterizing the range of the interaction between the 
side-loops in the effective potential centered on the back- 



bone monomers V gc to be w 



R g + 4cr = 6.95cr. This ef- 



fective interaction between the monomers was then added 
to the WCA potential, governing the non-bonded repul- 
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FIG. 5: (Color) Equilibrium structure for a circular backbone 
chain of length h = 400cr, with side-loops of length l s — 
20a attached to each backbone monomer. The polymer is 
confined within a cylinder of length L = 50.4<r and diameter 
D — 33.5(7. The splitting of the backbone into two parallel 
linearized branches is highlighted by using two color-codes, 
orange and blue, for the branches. The side-loop monomers 
are shown as transparent green beads. 



sive interactions, for a linear chain of length lb = 200a 
confined within a cylinder of diameter D — 29.5a and 
length L = 50.75<7. As FigfJJ} shows, this effective poten- 
tial reproduces the helical equilibrium structures of the 
polymer remarkably well. The structure factor displays 
a maximum at the same helical pitch value 4 /4 = 50a as 
the original simulations and reproduces the oscillations 
of the tangent-tangent correlation function (Fig. [2]) with 
only a slight global phase shift. The amplitude of the os- 
cillations (and hence that of the structure factor) is some- 
what larger for the effective potential, but this is prob- 
ably due to an overestimate of the interaction strength 
a. 

Finally we address what happens if the backbone is a 
ring-like polymer, appropriate to modeling, e.g., the cir- 
cular chromosome of E. coli. To that end we simulated 
a polymer with a circular backbone with Z& = 400<7 and 
side-loops of length l s = 20a, trapped within a cylinder 
of length L = 50 Aa and diameter D = 33.5c. The pack- 
ing fraction of the monomers is r\ = 18.9%, comparable to 
the one in the simulation of the linear backbone polymer. 
Fig. [5] shows that in this case the backbone loop is now 
organized into two parallel helices running along the long 
axis of the cylinder. As is evident from the snapshot, and 
corroborated by the analysis of the tangent-tangent cor- 
relations (Fig. [5] in Supplementary Material), the degree 
of helicity is reduced as compared to the linear backbone 
case due to the smaller side-loop length, but nevertheless 
remains significant. 

In conclusion, we have shown that the interplay be- 
tween the effective stiffness and intra-chain packing ef- 



fects caused by side-loops in polymers leads to novel 
helical equilibrium configurations of confined polymers. 
These structures are strikingly similar to ones recently 
observed in bacterial nucleoids. To what extent the phys- 
ical effects discussed here are sufficient to explain all 
the details of the large scale chromosome organisation 
in real bacteria is a question which clearly requires fur- 
ther research. At the very least, however, our results 
once again indicate that the ubiquitous aspecific interac- 
tions between the segments of long biopolymers like DNA 
can by themselves lead to significant spatial structuring, 
as has previously also been observed in the context of 
chromosome organisation in the nuclei of plants [24j and 
humans [25j . This argues for a more prominent place for 
polymer physics in the research into the structure and 
function of chromosomes. From a purely physical point 
of view our work points to novel possibilities for "sculpt- 
ing" the configurations of confined polymers by judicial 
choices of polymer topologies. 
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SUPPLEMENTARY: SPONTANEOUS HELICITY 
OF A POLYMER WITH SIDE-LOOPS 
CONFINED TO A CYLINDER 
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FIG. 6: The distribution of center of mass distance of the side- 
loops around the main-chain, n denotes the component of 
the distance along the local tangent of the main chain and r n 
denotes the component in a plane perpendicular to the local 
tangents. P(r n ) has a pronounced maximum near r n — 4a, 
whereas P(r t ) has a Gaussian shape with a variance of 8a 2 . 
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FIG. 7: The tangent-tangent correlation (inset) and its 
Fourier transform for a WLC polymer in presence of WCA 
repulsion between non-bonded monomers. 



We calculate the probability distribution of the center 
of mass position of the side-loops around the main chain 
polymer (Fig. This is done using our direct simu- 
lation of a main-chain polymer attached with side-loops 
confined within a cylinder. The main-chain has length 
lb = 200cr to each monomer of which is attached side- 
loops of length l s = 40(7. The polymer is confined within 
a cylinder of length L = 50.75(7 and diameter D = 29.5(7. 

The center of mass position of the loop r is measured 
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FIG. 8: The tangent-tangent correlation and its Fourier trans- 
form for a circular main-chain polymer attached with side- 
loops - the main-chain has length It = 200(7 and side-loops 
of length i s = 40(7 - and confined within a cylinder of length 
L = 50.4a and diameter D = 33.5c. The correlation (inset) 
shows clear periodic oscillation, and the periodicity is cap- 
tured by the peak in the Fourier transform at q p = 4. 



from the main-chain monomer to which the loop is at- 
tached. This position vector can be projected onto the 
direction of a local tangent on the main chain t, and 
a normal direction n such that r = r n n + r{t, where 
r n = r.n and r t — r.t. We calculate the distribution 
-P(rt), P(r n ) using 1000 equilibrated configurations of 
the polymer (Fig. EJ. 

We find that the center of mass of the side chains 
are distributed in a Gaussian manner in the direction 
parallel to the local tangent P(rt). Howeever, in the 
perpendicular planes, the distribution P(r n ) has a pro- 
nounced maximum at a distance r n = Aa. This off-center 
location of the center of mass of the side-loops, along 
with the radius of gyration R g = 2.95(7 leads us to use 
w ~ R g + 4(7 = 6.95(7 for the effective inter-side-loop 
interaction f3V gc {r) = a exp[— (r/w) 2 ]. 

In Fig. [7| we plot the tangent-tangent correlation and 
its Fourier transform for a self-avoiding (non-bonded 
monomers repel by WCA potential) WLC polymer of 
chain-length 4 = 200(7 and persistence length P = 
confined within a cylinder of length L = 50.75(7 and di- 
ameter D = 29.75(7 (corresponding to Fig. 4(a) of main 
text). Note the reduced magnitude and pitch of the helic- 
ity comapred to that of self-avoiding polymer in presence 
of the Gaussian-core repulsion as shown in Fig.l of main 
text. 

In Fig. [8] we show the tangent-tangent correlation and 
its Fourier transform for the case when the backbone is 
a ring-like polymer shown in Fig. 5 of main text. The 
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correlation function shows nice periodic oscillation which at q p — 4. This peak corresponds to a helical pitch of 
is captured by the maximum in the Fourier transform A ma:r = lb/A = 50 a. 



